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Abstract 

\£) • Existence of breather (spatially localized, time periodic, oscillatory) solutions of the topo- 

logical discrete sine-Gordon (TDSG) system, in the regime of weak coupling, is proved. The 
£NJ , novelty of this result is that, unlike the systems previously considered in studies of discrete 

breathers, the TDSG system does not decouple into independent oscillator units in the weak 
coupling limit. The results of a systematic numerical study of these breathers are presented, 
including breather initial profiles and a portrait of their domain of existence in the frequency- 
coupling parameter space. It is found that the breathers are uniformly qualitatively different 
from those found in conventional spatially discrete systems. 
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1 Introduction 



Soliton dynamics in spatially discrete nonlinear systems plays a crucial part in the modeling of 
many phenomena in condensed matter and biophysics. One also encounters such systems in the 
numerical simulation of soliton dynamics in high energy physics where they arise as approximations 
of continuum models. In this case one is concerned with the extent to which these discretizations 
resemble their continuum counterparts - of course one hopes for a close resemblance. On the other 
hand, for applications in condensed matter and biophysics it is crucial that in certain respects the 
discrete systems have features not found in their continuum counterparts. One aspect in which recent 
work has revealed such a distinction is in the appearance of so-called breather solutions - that is 
spatially localized, time periodic solutions. Such solutions are conventionally called breathers, by 
analogy with the breathers of the continuum sine-Gordon equation. While it is rare for a nonlinear 
partial differential equation to possess breather solutions (such solutions seem to be a feature of 
integrable systems) , there has been mounting evidence that they are ubiquitous in spatially discrete 
systems. After numerical work and some heuristic arguments for the existence of such solutions, 
the existence of breathers in a weak coupling limit of a class of discrete systems was proved by 
MacKay and Aubry [|j . Subsequently these ideas were extended to a wider class of discrete systems 
by Sepulchre and MacKay 0. The idea of both proofs is that each model has a parameter a, the 
coupling constant, and in the special case a — (called variously the anti-continuum, anti-integrable 
or weak coupling limit) the system decouples into independent oscillator units. If one sets one such 
unit oscillating while holding all others at rest, one obtains a trivial one-site (or unit) breather. 
Persistence of breathers in a sufficiently small neighbourhood of a = then follows from an implicit 
function theorem argument. 
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In this paper we present a study of breathers in the topological discrete sine-Gordon (TDSG) 
system ||. The system is so called because it reduces to the sine-Gordon equation in the contin- 
uum limit, and because it preserves many of the "topological" features of this equation, namely a 
topological lower bound on kink energy, and a continuous moduli space of static kinks saturating 
this bound. As a consequence of the preservation of these properties, kink dynamics in this system 
is remarkably close to that of its continuum counterpart, even in the highly discrete regime (kinks 
experience no Peierls-Nabarro barrier, and may propagate freely, indefinitely). One motivation for 
the present work is to see whether this continuum-like behaviour extends to breathers also. The 
TDSG system falls just outside the scope of the existence theorems actually proved in [jjj and |2j 
(although both papers indicate that their methods should apply in situations not directly covered 
by the theorems they state). The reason is that there is no limit in which the system decouples 
into independent units. Nevertheless, the analysis of the model and, in particular, its kink solutions 
suggest that there is an anti-continuum-like limit. For example, there are one-site breathers in this 
limit and one suspects that despite the lack of complete decoupling the continuation methods of 
MacKay et al. can be adapted to this setting. We show that this is indeed the case and prove a 
theorem about constant period continuation of the one-site breathers. 

Furthermore, the continuation can be performed numerically by searching for fixed points of the 
period mapping of the system using a Newton-Raphson method, a technique developed for other 
systems by Aubry and Marin Q . In this way we obtain a portrait of the period-coupling parameter 
space showing how far breathers of any particular period can be continued away from the zero 
coupling limit. We also obtain approximate breather initial profiles, and these display a rather 
unusual feature. Namely, the sites adjacent to the central oscillating site (the only one moving 
in the a = limit) shift away from equilibrium in the opposite direction to the central site. So 
the initial profiles are all, roughly speaking, sombrero shaped, rather than hump shaped as have 
usually been observed in discrete systems. The direction of continuation of breathers away from 
a = can be computed by a completely different (and very simple) numerical method, providing 
independent confirmation of this behaviour. The conclusion is, then, that breathers in the extremely 
discrete regime of the TDSG system are qualitatively quite different from those which appear in its 
continuum limit, and the two types of breather appear to be unrelated. 

The rest of this paper is arranged as follows. Section 2 introduces the TDSG system and derives 
the one-site breathers in its weak coupling limit. In section 3, persistence of these breathers away 
from zero coupling, and exponential spatial localization of the continued breathers are proved. Sec- 
tion 4 presents the results of the numerical search for breathers using the Newton-Raphson method. 
In section 5 the direction of continuation is rederived by independent means. The work is summarized 
in section 6. 



2 The topological discrete sine-Gordon system 

The TDSG system (introduced by Ward and first studied in || ) is a field theory in 1 + 1 dimensions 
where space is discrete, with lattice spacing h, and time is continuous. It is most conveniently defined 
in the Lagrangian formalism. To each two-sided sequence ip : Z —> R one associates the potential 
energy 

E P W = ^(Di + F*) (1) 

where 

2 1 1 

D n = t sin -(t/> n +i - ipn), F n = sin -(i> n +i +i>n)- ( 2 ) 

In the limit h — > 0, D — > ip x and F — > simp, so Ep reduces to 

1 r°° 

E P ° nt [i>] = - (^ + sin 2 V)dz, (3) 
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the potential energy functional of the continuum sine-Gordon model. The key feature of this par- 
ticular discretization is that 

D n F n = --(cOS"0n+l - COS^n), (4) 

so any configuration ip satisfying kink boundary conditions 

lim tp n = 0, lim ip n — n (5) 

n — > — oo n — >oo 

is subject to the inequality 

- \ J2( Dn ~ F «) 2 = Ep + \ Z]( cos ^n+i - cos^„) =E P -1. (6) 

So £"p > 1 with equality if and only if D n = F n for all n. Solutions of this first order difference 
equation for ip n are minimizers of Ep among all sequences satisfying the kink boundary conditions 
(||), and hence static solutions of the model. The general solution of the first order difference equation 
is 

ip n = 2tnn- 1 e a( - nh -^ (7) 

where a = h^ 1 log[(2+/i)/(2— h)] and b is a real valued parameter, identified as the kink position. So, 
rather remarkably, the TDSG system has a continuous moduli space of static kinks (parametrized 
by b £ R) and no Peierls-Nabarro barrier (Ep = 1 independent of b) . 

Dynamics is introduced into the system by defining the obvious kinetic energy 

Ek = \Y.&> (8) 

The action functional is then S[ip] = J (Ek — Ep)dt. The Euler-Lagrange equations associated with 
this functional, 

4 - h 2 4 + ft, 2 
= 4 , 2 cosV>n(smV>n+i + sin i/>„-i ) sini/>„(cos?/'n+i + cos-!/> n _i) (9) 

provide the equations of motion for the TDSG. The analogue of the anti-continuum limit of this 
discrete sytem is not h — oo, but rather h — 2, where the coefficient in front of the first term of (^|) 
vanishes, and where the static kink solutions (Q) degenerate to step-like functions. It is convenient 
to define a coupling constant 

4- h 2 

so that the "anti-continuum limit" is a = 0. In this limit the equations of motion reduce to 

■ipn = --sin^n(cos^ n+ i +COS^ n _i). (11) 

Note that the lattice does not decouple into separate oscillators, but maintains nearest-neighbour 
interactions. Nevertheless, as was observed in the a = system ([□]) has one-site breather 
solutions 

where 9(t) is any periodic solution of the pendulum equation 

(9 + sin (9 = 0. (13) 

We have chosen to locate the one-site breather at n — 0, but clearly any other site would do just 
as well. In fact, one could construct multisite breathers where any number of sites, none adjacent 
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to any other, oscillate according to the pendulum equation, with periods in rational ratio. In this 
paper we will consider only one-site breathers, however. 

As a system of autonomous differential equations, there is an IR-action on solutions of the TDSG 
system corresponding to the time translation symmetry of the system. Thus any breather solution 
gives rise to a 1-parameter family of breathers. As a consequence, a breather solution is never an 
isolated solution J5j. To apply an implicit function theorem argument one needs to deal with this 
degeneracy - on the full phase space of the system the linearization of the equations of motion at 
a breather always has a non-trivial periodic solution which generates time translation along this 
breather. We deal with this by considering a subset of the phase space corresponding to solutions 
with a certain time reversal symmetry. We choose to work with solutions which are odd with respect 
to the time reversal symmetry i.e. solutions that satisfy t) = — cf. |Q]. The main reason for 
this choice of symmetry is simply that odd solutions of the pendulum equation are simpler to write 
down in terms of elliptic functions than are even solutions. We will then prove that the continuation 
theorem in the space of odd functions implies a similar result in the space of even functions. Were we 
to consider continuation of multisite breathers, we would need to remove several phase degeneracies 
(one for each excited site) in order to apply the implicit function theorem. In fact the same procedure 
(restriction to odd solutions) removes all such degeneracies. 



3 Continuation of one-site breathers: analytic results 

Our proof of the existence of discrete breathers for the TDSG system in the weak coupling regime 
follows the method used by Aubry and MacKay in [Q. The proof proceeds in two steps; the first 
shows that certain periodic solutions persist and the second that these solutions continue to be 
exponentially localized. In this section, after introducing the function spaces we choose to work 
with, we state and prove the two halves of the continuation theorem. 

3.1 The function spaces 

Let C^ 2 denote the space of C 2 real- valued T-periodic functions on R with odd time reversal 
symmetry. Equipped with the standard C 2 norm 

Ha = sup max {^(i) |, |^(t)|, |^(i)|}, (14) 
ten 

C^ 2 becomes a Banach space. Similarly, let C^ denote the space of C° real- valued T-periodic 
functions on R with odd time reversal symmetry, equipped with the standard C° norm 

Ho=sup{|V(t)|}. (15) 

teR 

Consider the linear space ©„ e zA„, where each X n = C<^ 2 , equipped with the natural sup norm i.e. 

if * = (^n)nez € (SnezX n then 

||*|| 2 =sup|* n | 2 . (16) 
Definition 1 Qy 2 := {* e 0„ eZ A„ : ||*|| 2 < oo} 

Similarly, equip the linear space ® n eiXn, where each Y n = C^ , with its natural sup norm i.e. if 

* = (VvOnez e ®neiXn then 

||*|| = 8up|* n | . (17) 
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Definition 2 Q TQ := {* e ® n &Xn ■ ||*||o < 00} 

Both f2^ 2 and are Banach spaces. 

The equations of motion (^|) give a natural map F : 2 © ^ — * defined by 

F(9,a) n = i>n - acosVVi(sinV>n+i +sin^n-i) + (a + -) sin<0n(cosV>n+i +cos^„_i). (18) 

The zeros of F correspond precisely to the odd T-periodic solutions of the equations of motion (0) . 

In section |^ we argued that a — should be considered as something akin to an anti-continuous 
limit for the TDSG system. Although the equations of motion ([□]) do not (in the terminology of 
MacKay and Sepulchre) decouple into local units, the energy density of the kink solutions becomes 
concentrated at one site - as one would demand for an anti-continuous limit. Moreover, like the 
models considered by Aubry and MacKay, at a = the system admits single site breather solutions - 
in our case the single non-stationary site moves according to the pendulum equation. So, analogously 
to the continuation of one-site breathers from the anti-continuous limit considered by various authors 
[fjj ||, one can ask what happens to these solutions as the coupling a is increased. Despite the 
somewhat different nature of the TDSG system we find similar continuation results to those of the 
authors mentioned above. 



3.2 Statement of the continuation theorem 

Given T S (2tt, 00) \ 2nN, denote by 9 the unique odd T-periodic solution of the pendulum equation 
( p"3| ) with positive initial velocity. Let denote the corresponding one-site breather solution of the 
TDSG at a = given by JT^). Then we have the following theorem. 

Theorem 3 There exists e > such that for a € (0, e) there is a unique continuous family, ® a , of 
odd T-periodic solutions of TDSG at coupling a, such that &q = 0. Moreover, these solutions are 
exponentially localized, i.e. there exist positive constants B and C (depending on T and a) such that 
|(0a)nlo ^ C exp(— 5|n|) holds for all n G Z. 

In other words, for small enough a the one-site breather © has a locally unique continuation among 
odd T-periodic solutions and these continued solutions remain exponentially localized. The proof 



of this theorem is given in the next two sections. Section 3.3 establishes the continuation through 



T-periodic solutions while section 3.4 deals with the question of localization. 



3.3 Persistence of periodicity 

By our remark in section [O], is an odd solution of the TDSG system with period T for coupling a, 
if and only if, F(\P,a) = for the function F defined by equations (|l8|). The odd one-site breather 
of period T gives us a solution for a — 0. The proof proceeds by showing that F is a C 1 function 
to which the implicit function theorem may be applied. Thus a locally unique zero (and hence an 
odd T-periodic solution of TDSG) of F for a close enough to is demonstrated. It remains to show 
that the implicit function theorem can be applied to F, i.e. we need to check that F is C 1 and that 
DF, its derivative with respect to the fi^ 2 factor is invertible at the one-site breather solution. 

Proof: The verification of the first fact is straightforward. A short computation gives us the following. 
If we let 

DF : (Q~ 2 © R) (0" 2 © R) — ► fl~ (19) 
be the derivative of F and write 

DT*, tt (,5*,fo) = (^ n ) nez (20) 
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then DT is defined by 

tp n = dipn — [sin^n (cos^„+i +cosip n -i) — cos<0„ (sin^ n+ i +sin^„_i)]<kt! 

+ [asinipn (sinV^i+i + sin^>„_i) + /3cos-0„ (cosV'n+i + cosV'n-i)] 5ip n 
-acos?/; n (cos ip n +i Stp n +i + cos-0„-i 8tp n -i) 

+/3 sinip n (smip n+1 Sip n+1 + sin^„_i 6ip n -i) , (21) 

where (3 = a + h. From this it is easily seen that DT is bounded. 

At a = 0, £)F, the derivative of F with respect to the factor, reduces to 

DF^(S^) n = 8ti>„ + ^ cosi/) n (cos V'n+i + cos tp n -i)5ip n 

+ i sinV>™(smi/>n+i<5V> n +i + sin ip n -±5ip n -i) . (22) 

In general this is a tridiagonal but not diagonal operator. However, DF evaluated at any one-site 
breather is diagonal. Explicitly, substituting VP = 0, the one-site breather of the form (jT^), into 
(|2|) yields DF e = 0„ eZ T n , where L n : C^ 2 -> C^g is, 

Jf/) + cos 9ip if n = 0, 

Vi + i(l + cos6»)^ ifn = ±l, (23) 
i/j + ip otherwise. 

Clearly, to prove invertibility of DF at a one-site breather it is sufficient to prove that each operator 
L n is invertible, since then DF^ 1 — (BnezL^ 1 and we see t hat H- DT -1 ! < max{|TQ |, \, IL^ 1 ]}- 



Invertibility of the operator L n will follow (see section 3.3.2 ) from standard facts about linear 
differential equations once we establish that kerL„ = {0}. 

3.3.1 Injectivity of DF e 

For T £ 27rZ, one sees immediately that kerT„ = {0} for \n\ > 1. So it remains to consider the 
central (n = 0) and off-central (n — ±1) equations 

+ = (24) 

-(l + cos0)V> = (25) 

where 9 is a T-periodic, odd solution of the pendulum equation, 9 — — sin#, with positive initial 
velocity. Both (|2~f| ) and (|25| ) are Hill equations, that is ODE's of the form %i> = q(t)ip where q 
is periodic. Such equations have a long history. They arise, for example, in stability analyses of 
celestial mechanics M, Note that in both (^J) and (|2f]) the potential, q(t), has basic (i.e. smallest) 
period T/2 since cos is even. Such equations may have periodic solutions only with basic period not 
less than T/2. To prove that Lq and L±\ have trivial kernel, we will have to prove separately the 
non-existence of T/2 and T-periodic odd solutions of ( pij ) and (|25|). 

This task is simplified slightly for equation ( p4[ ) by the observation that ip = 9 is an even T- 
periodic solution. This is, in fact, a consequence of the time translation symmetry of one-site breather 
solutions previously mentioned. A standard result of Floquet theory Q implies that ([m]) cannot 
have both T/2-periodic and T-periodic solutions, so it suffices in this case to rule out existence of 
T-periodic odd solutions. No such simplification occurs for the off-central equation (|2^). 

To proceed further one must consider the specific form of the solution 9 of the pendulum equation. 
If 9(0) — and 9(0) — 2k, where k G (0, 1) then the solution has amplitude of oscillation 9o G (0,7r) 
satisfying 9q = 2arcsinfc, and period T = AK where 

K = r , ^ (26) 
Jo Vl-& 2 sin 2 z 



G 



is the complete elliptic integral of the first kind ||. The parameter k is called the modulus in the 
literature on elliptic functions. In fact this solution satisfies 



sin = ksnt 



(27) 



where sn is the sn-oidal elliptic function of Jacobi. Note that sn itself, like all elliptic functions, de- 
pends parametrically on k. This dependence will always be suppressed in our notation. Substitution 
of ( |27| ) into ( |2~i| , |2"5| ) recasts them as a pair of Lame equations (9) , 



^+ (1 - 2k 2 sn 2 t)tp = 
ij> + (1 - k 2 sn 2 t)ip = 0. 



(28) 
(29) 



A standard technique when seeking 2K and AK -periodic solutions of a Lame equation is to make 
the following elliptic change of time coordinate, 



1 1 * x(t) = am(i) 



(30) 



where am is Jacobi's amplitude function, that is, the unique continuous function satisfying sinx(t) = 
sni, x(0) = 0. Under this coordinate change the pair of Lame equations above becomes a pair of 
Ince equations [|], 



(1 + acos2x)ip" — as'm2xip' + (I — a + 2acos2x)ip = 
(1 + acos2x)ip" — asm2xtp' + (1 + acos2x)ip = 



(31) 
(32) 



where ip' denotes dip/dx etc., and a = k 2 /(2 — k 2 ) is a period dependent parameter taking values 
in (0, 1). Note that the coefficients in these ODE's are periodic in x with period n, independent of 
k, and hence T. So all k dependence is now explicit and resides in the constant a(fc). Odd 2K and 
4i^-periodic solutions of (28|29|) correspond to odd 7r and 27r-periodic solutions of (|3l]J3l), and vice 



So, we need to prove that ( |3ll) has no odd 27r-periodic solutions, and that (^2|) has neither 2ir 
nor 7r-periodic odd solutions. Assume, to the contrary, that such solutions exist. Let ip , ip 2 be 
27r-periodic odd solutions of (|3l]) and (|3^ ) respectively, and ip 3 be a 7r-periodic odd solution of 
([52J) . Functions of such parity and periodicity must have purely sinusoidal Fourier expansions of the 
following forms: 



v> 3 



^2 K. sin( 2n - £7 = 1,2 

n=l 
00 

^2 A n sin 2nx - 



(33) 



n=l 



Since the coefficients in ( |31| ) and (32) are analytic in x, each ip a (x) is analytic in a strip of the 
complex x plane containing the real axis |l(|. It follows that the Fourier coefficients A°, a — 1,2, 3, 
must satisfy strong decay criteria, namely 



lim n v A a n = 

n — >oo 



(34) 



for every p > 0. 

Substituting the Fourier expansions for ip a into their appropriate Ince equations one obtains the 
following linear second order difference equations for A^: 



A\+A\ 



0, 
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A 1 



An 2 — An + a 



A 3 



o(2n-l)(n + l) " 
8n(n- 1) 



n-2 
n — 1 



- A„ 



A? + 5A 2 2 
-A„, 



a(An 2 + 2n- 1)"" ' 4n 2 + 2n - 1 " n ~ l 

6A? + WaA\ 
An 2 - 6n + 1 



2(4n 2 - 1) 
a(4n 2 + 6n + 1) 



4n 2 + 6n + 1 



4 3 

A n-1 



n > 2; 
0, 

n > 2; 
0, 

n > 2. 



(35) 



(36) 



(37) 



Lemma 4 For a/Z a G (0, 1), aZ/ nontrivial solutions of the difference equations ( p^ , pq . |^ j are 
exponentially divergent. (By trivial we mean identically zero.) 

The proof of this lemma is straightforward, but rather involved. The interested reader is directed 
to the appendix. 

Since none of the sequences A* satisfy the decay criterion (|34|), we conclude that no such solutions 
ip a exist. Hence kerLg = kerL±i = {0}. 



3.3.2 Invertibility 

Having established that each L n as a linear operator C^ 2 — > C^ is injective, we apply standard 
linear ODE theory to show that each L n is invertible (i.e. its inverse is bounded). Consider L n 
acting on the full space Ct.2 — ► Cr,o- Each L n respects the splittings Ct,* — C^ * © Ct*> where 
Cj, # denotes the subspace of even loops. 

Lemma 5 keri„ C C^ 2 . 

Proof: This is immediate from previous work except for the cases \n\ — 1. Suppose ip = ip + + 1\>" G 
Ct.2 and L n ip = 0. Then since L n preserves the odd-even splitting, we must have L n tp + = L n ip~ = 
and hence ip~ = 0, by the previous section. □ 

Let keri n = {ij; G Ct.o '■ (Vs kerL n ) 2 = 0}, where ( , ) 2 is the standard L 2 inner product. 

Corollary 6 C^ Q C keri n . 



Proposition 7 Each L n : C T 2 — > C T is an invertible linear map. 

Proof: Standard linear ODE theory implies that L n : Ct,2 —> kerLn 1 " is surjective. In particular, 
L n : Ct.2 - > Ct o kerZ n is surjective. But since L n preserves the odd-even splitting, L n : C^ 2 — > 
is still surjective. Now we have a bounded, injective, surjective linear map and hence the closed 
graph theorem implies that L n has a bounded inverse. □ 

3.4 Exponential localization 

We now turn to the question of localization. Issues related to exponential localization in the general 
class of networks have been treated by Baesens and MacKay in [jll| . Localization of our continued 
periodic solutions will follow from an application of their results. To state the relevant results from 
[ p"i"| we need to recall some definitions from that paper. 

Consider a network, that is a countable metric space S with metric d (in our case S = Z with 
d(m,n) = \m — n\). To each s e S associate two normed spaces X s and Y s , with norms |.|x, and 
|.|y s , to be considered as a local state space and a local force space respectively - in our case each 



8 



X s is the space of C 2 odd T-periodic loops in K with the norm given by equation (|l4|) and each Y s 
is the space of C° odd T-periodic loops in R with the norm given by equation (|lj]). Define X by 



X = G (B S X S : \\x\\ x := sup |a; s | Xa < oo j (38) 

and Y similarly (in our case X = 2 , Y = 17^ and the norms on X and Y are just those already 
given on these spaces of loops in section 3). 

Definition 8 A bounded linear map L : X — > Y is said to be ^-exponentially local if there exists 
e > and a continuous function <p : [0, e) — ► M. such that 

sup V \L rs \z d ^ < 4>[z), z G [0, e). (39) 

The tridiagonal linear map DFq : 2 — > is an exponentially local map (e.g. we can take 
<f>(z) = \\DF&\\(2z + l)). 

Definition 9 u G Y (or in X ) is (C, A)-exponentially localized about a site o G S, if there exist 
C > 0, < A < 1 such that \u s \ < CX d ^ s '°\ 

Certainly, a one-site breather is exponentially localized about the single oscillating site o (e.g. we 
can take C — \u Q \, and any A > 0). We now state the main result of fll|| . 

Theorem 10 (Baesens-MacKay) If L : X — > Y is <f>- exponentially local and invertible, and u G Y 
is (C, X)- exponentially localized about site o G S, then L~ l u is (WC, n)- exponentially localized about 

0, where W and fi depend only on ||L _1 |j ; A and a function (3 related to the function 4>. 

As an application of the above theorem, MacKay and Baesens prove a theorem concerning 
continuation problems. Let F : M. (B X ^ Y be a. C 1 function, which at a — has a regular zero, 

1. e. there is a solution Qq of F(Q, 0) = at which DFq is invertible. Then, by the implicit function 
theorem there is a neighbourhood of in which F(a, O) = has a unique solution Q close to Qq. 
The solution satisfies 

de a T -x(dF \ 

L a 37" (® a ,a) (40) 



da a \da 

where L a — DF a (Q a ) and may be continued as long as L a remains invertible. In this setting they 
prove: 

Theorem 11 (Baesens-MacKay) Suppose Qq is exponentially localized about site o G S, with con- 
tinuation Q for a G [0, ao] and with L a cf>- exponentially local. If 4^(Bo) is (A, fi)- exponentially lo- 
calized about o, F G C 2 , and -j^DF is if) -exponentially local for |0| < sup Q6 j 0QO ] \Q a \ and a G [0, ao], 
then ©q, is exponentially localized about o for a G [0, ao] • 

We have shown that for the function F defined by (|is|), for each T G (2tt, oo) \ 27rN the odd 
one-site breather, O of period T is a regular zero of F. Choose ao > so that for a G [0, ao] the 
continuation & a exists. We claim that the conditions of Theorem [ll] obtain. 

For each a G [0, ao], L a is ^-exponentially local where <fi(z) — A(2z+1) and A = max Qe r 0)Qo i ||iv a ||. 
From 

,„ f 2 sin6> if n = 0, 

^(©o)n = < -sin0 if \n\ = 1, (-11 i 

'" otherwise 



we see that ^(0o) is exponentially localized about n — 0. It follows from 

^-DF^(S^) n = [(COS (VWl - ^n) + COS(^„_l - tpn)) Hn 

da 

+ COS (lp n+ i + lp n ) Stpn+i + COS (lp n -l + i> n ) Hn-l] (42) 
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that -^DF is ^-exponentially local for ip(z) = 2(1 + z). It is straightforward to check that F is C 2 . 
Hence Theorem |ll| may be applied in our situation: the continued periodic solutions Q remain 
exponentially localized and thus are discrete breathers as claimed. □ 

3.5 Continuation of even breathers 

We have proved that for any T £ (2tt, oo)\27tN there exists a locally unique continuation ® a through 
fiy 2 of the odd one-site breather 0. The purpose of this section is to prove that a similar result 
holds for continuation of even one-site breathers also. Since this is not quite so obvious as it may 
seem, we will state and prove this result carefully. Let fir,™ be the space of all T-periodic C n loops, 



and be its even subspace, both equipped with their natural norms (see section 3.1). Translating 
0, the odd one-site breather of period T, through a time of T/4 one clearly obtains an even one- 
site breather, call it (so explicitly, 0(i) = &(t + T/4). Note this is a slight abuse of notation: 
we are using to represent both a point in 0^ 2 and the value of the corresponding function at 

a particular time). If one similarly translates the family Q , one obtains a continuation a of 
through breather solutions in Ot,2- We will prove that, for a small enough, this continuation 
actually lies in the even subspace 2 . 

Corollary 12 For each T £ (27r, oo)\27rN there exists e > such that there is a locally unique 
continuous family & a £ ilj, 2 of solutions of the TDSG system at coupling a £ [0, e), with ©o = 0. 
These solutions are exponentially localized. 

Proof: By Theorem |^ there exists e > such that for a £ [0, e) there is a continuous family 
©a £ 57^2 °f °dd breather solutions with ©o = ©. Choosing e sufficiently small, each Q is 
unique (i.e. is the only zero of F) in the open ball B e (&) C U T2 . Consider the family @ Q defined 

by & a (t) — ® a (t + T/4). Clearly Q is a continuation of the even one-site breather © through 
exponentially localized breathers in Qt,2- It remains to prove that this continuation actually lies in 

For each a £ (0, e) define another breather solution & a of the TDSG system by & a (t) = 
—® a {t + T/2). By periodicity and oddness of @ Q one sees that & a £ f^ 2 - Also, — @|| 2 = 
||0 Q — @|| 2 < e since the one-site breather satisfies the identity 0(i + T/2) = — &(t). But 
the continuation Q is unique in the ball _B e (0) c f^ 2 ' so we conclude that Q = ® a . Hence 
& a (T/2) = -00,(0), and by time reversal symmetry of the TDSG system, a (T/4) = 0. But 

Q (O) = © a (T/4) = 0, so Q is even. 

Local uniqueness of the even continuation Q is immediate. For if the continuation were not 
unique on a small enough neighbourhood of 0, one would be able to construct an alternative odd 
continuation, again by time translation, in violation of local uniqueness of Q . □ 



4 Continuation of one-site breathers: numerical results 

For each non-resonant period T, Corollary [l2] yields a 1-parameter family of T-periodic even breathers, 
but fails to give detailed information about these solutions. For example, it gives little information 
about the domain of existence of the continuation, how this domain varies with the period T, or the 
profiles of the breathers. We have obtained some such information through numerical work. This 
section describes in detail the numerical scheme used and the results obtained. 

Before describing the numerical work we note that there is one established mechanism which puts 
limits on the extent of the continuation being attempted. Namely, one expects unique continuation 
to fail once the breather enters into resonance with the phonons of the system. Linearizing the 
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Figure 1: The continuation domain of one-site breathers in (w,a) space. The upper, jagged curve 
is the boundary obtained from the phonon resonance argument, and consists of parabolic segments. 
The shaded region beneath this curve is the actual continuation domain obtained by numerical 
analysis for ui E [0.25,1], one vertical line for each frequency sampled. The region to the left of 
oj = 0.25 is inaccessible to our numerical scheme. 



equations of motion about the vacuum and seeking solutions of travelling wave type yields the 
phonon dispersion relation 

k 

uj 2 = 1 + 4a sin 2 - (43) 

relating frequency to to wavenumber k ||. From this we see that the phonon frequency band at 
coupling a is [1, \/l + 4a]. Suppose we fix a frequency u> £ (n^ 1 , (n — 1) _1 ), n = 2, 3, . . ., so that 
its lowest harmonic above the bottom edge of the phonon band is nuj. Then during continuation 
of the frequency u> one-site breather, a increases from and the phonon band expands, eventually 
capturing the harmonic ntu when nuj = vl + 4a (if continuation persists this far). Beyond this point 
one expects not only loss of uniqueness of the continuation, but also a breakdown in exponential 
localization (l2), since the phonons themselves are not localized. Thus the union of the curves 

a = \{n 2 u? - 1), rT 1 < oj < (n - l) -1 , n = 2,3, ... (44) 

with vertical line segments at the resonant frequencies ui — rT 1 provides an upper boundary to the 
possible breather existence domain in the (to, a) parameter space (see figure 1). We will find that 
the existence domain is considerably smaller than the region bounded by this curve. 

The idea behind the numerical method used to map out the existence domain is to convert the 
existence proof of section ^ into a constructive means of finding periodic solutions. By the analysis 
of that section, we know that the nonlinear function F defined therein is invertible at a one-site 
breather (provided the frequency of the breather is nonresonant) and hence, by the implicit function 
theorem, we know that in a neighbourhood of the one-site breather, _F -1 (0) is a smooth curve. 
This curve is the one-parameter family of periodic solutions we seek, parametrized by the coupling 
a. Given one point on the curve (e.g. the one-site breather at a = 0) one can find sufficiently 
close neighbouring points by using the given point as initial input in a root-finding algorithm for F. 



11 



Working piecemeal from a — 0, we can in this way construct the entire curve. To solve F — we 
may use a Newton-Raphson method, which will converge provided DF is invertible at the root and 
we attempt a small enough step along the curve. 

As it stands, this scheme is impractical since the spaces on which F is defined are infinite 
dimensional. We may seek a finite-dimensional truncation of the scheme in several ways. One 
approach is to truncate the lattice to finite size, and approximate the space of periodic functions 
associated with each lattice site by some finite Fourier series representation. Then the problem is 
reduced to that of determining the finite number of Fourier coefficients at the finite number of lattice 
sites, which is solvable by Newton-Raphson. This is somewhat complicated in the TDSG system 
because one must compute Fourier transforms of non-polynomial nonlinear terms, necessitating the 
use of fast Fourier transforms [[[3] . 

For this reason we prefer an alternative approach, where F is replaced by a function stemming 
from the Poincare return map for period T, which we shall denote T. We again truncate the lattice 
to finite size (say 2 N + 1 sites symmetrically distributed about the breather centre) , and consider 
the map on phase space which maps the initial data of a solution to its phase space position after 
time T. This map may be computed in practice by solving the lattice equations of motion using a 
4th order Runge-Kutta method. Fixed points of the map T, or equivalently, zeros of T — Id, are 
clearly periodic solutions. However, on the whole phase space, solutions are never isolated because of 
the time translation symmetry. Consequently, DT is never invertible, and hence we cannot directly 
apply a Newton-Raphson algorithm. 

The solution to this difficulty is to consider some restriction of T which has no such time- 
translation symmetry. Marin and Aubry Q consider the restricted return map Tp : R 2Ar+1 — > M 2Ar+1 
which takes initial positions of time-reversal symmetric solutions (so that initial velocity is 0) to their 
positions at time T. That fixed points of Tp are truly T-periodic is an immediate consequence of 
energy conservation. An unfortunate consequence of this choice is that at a = 0, D(Tp — Id) is 
singular for all lu. Explicitly, 

-DTp = diag(cos T, . . . , cos T, Vl (T) , y (T) , Vl (T) , cos T, . . . , cos T) , (45 ) 

where yo,yi are solutions of ( |24|j25| ) respectively with initial data y n (0) = l,y n (0) = 0, 6 now 
being the even T-periodic solution of the pendulum equation. Let yo be the odd solution of ( |24| ) 
with y a (0) = 1. Then a short calculation reveals that yo(t) = —0(t)/ sin 6(0) which is T-periodic. 
Conservation of the Wronskian yoy — yoyo then implies that yo(T) = 1. Hence det(DTp — I) = 0, 
and one is not guaranteed good behaviour of the Newton-Raphson algorithm for Tp — Id at, and 
close to, a = 0. This problem is not specific to the TDSG system, and was in fact pointed out in 
I- 

An alternative restriction, which we denote Ty, arises by mapping the initial position of a time- 
reversal symmetric solution to its velocity after one period. Certainly, T-periodic solutions are zeros 
of Ty , although the converse is clearly false since the 2T-periodic one-site breather is a zero of Ty at 
a = 0. (In fact any zero of Ty must be 2T-periodic by time-reversal symmetry of the system.) This 
is not a problem provided the zero we seek is isolated, which is guaranteed by the inverse function 
theorem so long as -DTy remains invertible along the continuation curve. The advantage of this 
approach is that at a = 0, 

DT V = diag(- sin T, . . . , — sin T, y x (T) , y (T) , m (T) , - sin T, . . . , - sin T) , (46) 

which turns out to be nonsingular except at the resonant periods T = 2?i7r, where continuation is 
not expected anyway, and also at the "false resonances" T — (2n + l)ir. As continuation proceeds, 
care must be taken, once a zero of det DTy has been encountered, that we do not follow a spurious, 
non-periodic branch of T^ 1 (0). At the end of the Newton-Raphson algorithm we check separately 
that the proposed periodic solution really is a fixed point of Tp to within numerical tolerance (note 
that this involves no extra computational cost since Tp is an inevitable by-product of the calculation 
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Figure 2: Details of the continuation domain around three resonant frequencies: (a) u) = 1/4, (b) 
u) = 1/3 and (c) lu — 1/2. In each case, the solid curve with crosses represents the numerical results, 
the dashed curve is the phonon resonance boundary for an infinite lattice, and the six solid curves 
without crosses are the parabolic resonance curves for the 13-site lattice. 

of TV). So it is possible that the continuation could fail prematurely at a zero of DTy if it is diverted 
from the true curve onto a spurious branch. This turns out to be a problem only around the false 
resonant frequency ui = § in our work (see discussion below). 

To test whether the continuation has failed we use two criteria, beyond mere failure of the 
Newton-Raphson algorithm to converge. First, we check that each new solution does not differ 
greatly from the previous one, by computing ||0 Ql — ® Qo ||2- This is to ensure that the continuation 
has not jumped across to some distant zero of Ty, as Newton-Raphson output is wont to do close to 
singular points. Second, we seek bifurcations of the continuation by applying a simple, necessary, but 
not sufficient, test. Namely, we terminate continuation if ever det DTy and dctDTp simultaneously 
vanish (note that this test requires minimal extra computational effort since DTy is required by 
the Newton-Raphson algorithm, and cannot be calculated without producing DTp as a by-product. 
The only extra cost is that of the LU decomposition of DTp, insignificant compared to the cost 
of the Runge-Kutta method). That this test is not sufficient is clear from consideration of the 
false resonant frequencies u> = — at a — 0, where no bifurcation can exist by Corollary A full 
bifurcation analysis, along the lines of Jl3| | , would remedy this insufficiency, but only at considerable 
computational cost, and would be well beyond the scope of this paper. We hope to perform such an 
analysis in the future. 

The results of the continuation algorithm, using a 13-site lattice with fixed endpoints and the 
implementations of the Runge-Kutta and Newton-Raphson methods outlined in |l4|], are presented 
in figure 1. The phonon resonance boundary clearly curtails continuation just to the right of the res- 
onant frequencies u> = -g, | , However for large regions of a; space, continuation fails before phonon 
resonance can occur, at what appear to be genuine bifurcations (at least detDTy = detDTp = 0) 
the underlying physical mechanism for which remains mysterious. The exception is a small u> range 
around the false resonance oj = | where the continuation seems to fail spuriously due to its en- 
countering a zero of det DTy which is not a genuine bifurcation point (detDTp ^ 0), but which 
derails the Newton-Raphson algorithm nonetheless. The u> € [0.27, 0.3] part of figure 1 should thus 
be viewed with a little scepticism. 

Closer examination of the results around uo — \t\,\ (figure 2) reveals that the continuation 
slightly oversteps the resonance boundary defined by equation (Q) . This can be well understood as 
a finite size effect: equation (^) is the dispersion relation for an infinite lattice, whereas we have a 
27V + 1 site lattice, divided by the central site into two "tails" of length N which remain close to the 
vacuum. The relevant spectrum for small amplitude oscillations is doubly degenerate, consisting of 
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Figure 3: Initial profiles of breathers, (a) at frequency 0.6 with coupling h — 2.0 (x), h = 1.8 (o) and 
h = 1.73 (+); and (b) at coupling h — 1.73 with frequency u = 0.8 (x), tu = 0.52 (o) and oj = 0.44 

(+)• 



two copies of the normal mode spectrum of the iV-site linearized lattice (one copy for each of the 
two tails), which itself consists of N distinct frequencies ujj, j — 1,2, ... ,7V, easily shown Q to 
satisfy 

uj] = 1 + 4a sin 2 (— j —Y (47) 

So bifurcation should be expected at that value of a for which ujn coincides with the lowest harmonic 
of uj above 1. This leads to a slightly different boundary from that previously discussed. Since the 
"phonon" spectrum is now discrete, it is possible for the continuation to skip over this first bifurcation 
point occasionally, and detect a later one instead, due to u>n-i, etc. In our case, N = 6, so 

there are six parabolic resonance curves, one for each normal mode frequency. These are included 
in figure 2. 

As a by-product of the numerical work described above we also obtain breather initial profiles 
within the continuation domain. The resulting profiles are highly spatially localized, and remarkably 
uniform, all exhibiting basically the same shape: compared to the one-site breather of the same 
period, the central site of the continued breathers is displaced slightly more than in the one-site 
breather, the two adjacent sites are displaced negatively by a small amount, and the remaining 
sites are barely displaced at all. Figure 3 shows breather profiles of the same frequency at increasing 
coupling, and profiles at the same coupling and various frequencies. Breather shape will be discussed 
further in the next section. 

Finally, we remark that the numerical methods described above are inefficient for investigating 
the long period breathers of the system. Indeed, if one used the same time step in the numerical 
integration schemes independent of the period, T, then of course the time to compute Ty and DTy, 
which occupies the vast majority of the computing time, would grow linearly with T. If one wanted to 
investigate long period breathers then it would seem more efficient to switch from the Runge-Kutta 
scheme to a symplectic integrator, of which one expects reliable long-time integration at modest 
time step size. The efficiency of the algorithm could also be improved by enabling it to take larger 
steps in a so that the bifurcation edge is reached more quickly. One method to do this (resulting in 
a so-called Euler-Newton scheme) is to generate the initial guess for the Newton-Rapshon algorithm 
by extrapolating along the tangent vector to the continuation curve at the last established point. 
Calculating this vector involves solving an extra set of coupled linear ODEs, but it is thought that 
the benefit of an enlarged parameter step length generally outweighs this extra computational cost 
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5 The direction of continuation: a linear calculation 



Corollary 12 shows that for each even one-site breather of nonresonant period T there exists a 
family & a € fi^ 2 , for a in a neighbourhood of 0, satisfying the continuation equation 



F(@ a ,a) = 0, 



(48) 



with ©o = 0. In this section we describe a simple method for computing the direction of continuation 
away from 0, that is d& a /da\ a= o € £^ 2 j which will henceforth be denoted 0' = (xn)nez- The 
method has previously been used to provide the "predictor" step in predictor-corrector schemes 



although we shall require only a greatly simplified 
Differentiating equation ( pl8| ) with respect to a at 



DF&&' = 0. 



(49) 



for continuation problems mentioned above [|16| 
version since we consider only the case a = 0. 
a = yields 

dF 

Denoting the T-periodic even solution of the pendulum equation 9, one easily obtains the following 
explicit expressions, 

2 sin 9 n = 

-sin0 n = ±l (50) 
\n\ > 1, 



OF 



and 



[DF e ®'] 




COS 9Xn 
±(l+COS0)x> 

Xn 



n - 
n ■ 
\n\ 



= 
= ±1 
> 1. 



(51) 



So inverting DFq to find 0' amounts to solving an infinite set of ordinary differential equations, 



Xo + cos 6»xo 

X±i + oC 1 + cos0 )x±i 



= -2sin6> 
= sin 6* 

= M > 1, 



(52) 
(53) 
(54) 



in the space f2^ 2 - Provided T ^ 2-7rZ, equation (]54|) has only the trivial solution \n — 0, so to first 
order the continuation leaves the \n\ > 1 sites fixed at 0. 

For the < 1 equations ([52],^3]), the question remains how one should choose initial data 
Xn(0) € K, Xn(0) = in order to generate a T-periodic solution. Let y n (t) be the solution of the 
homogeneous problem with initial data y n (0) = 1, y n (0) = 0, and Y n (t) be the particular integral of 
the inhomogeneous equation with initial data 1^(0) = Y n (0) = 0. Then 



Xn(t) = Xn(0)yn(t)+Y n (t) 

and demanding that Xn be periodic gives an overdetermined pair of linear equations, 

Xn(0) = y n (T) Xn (0) + Y n (T) 
= y n (T) X n(0) + Y n (T) 



(55) 



(56) 
(57) 



for x«(0) m terms of the constants y n (T) 1 y n (T),Y n (T),Y n (T). Equation ( p7\ ) can be solved for 
Xn(0), leaving ( p6| ) as a linear consistency condition for the constants. Clearly X-i(0) = Xi(0)- 

So the direction of change of the initial profile of a breather can be deduced by solving a coupled 
set of five second order ODE's (the pendulum equation, and the homogeneous and inhomogeneous 
equations for n = and n = 1), over the time interval [0, T], computing the end-point constants 
and solving (|57]). This has been done numerically using the fourth order Runge-Kutta scheme im- 
plemented by Maple. The results are presented graphically in figure 4, a plot of Xo(0) and x±i(0) 



15 



0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

to 

Figure 4: The components xo and \i 01 the tangent vector 0' as functions of frequency, u>. The 
solid curves were generated using the linear method, the dashed curves from the breather profiles 
generated by the numerical scheme described in section |4|. The two data sets agree closely (especially 
for xi). 

against frequency ui = 2tt/T. Note that for all u>, \o > and \±i < ; so the continuation starts 
by pulling the two off-central sites downwards, and the central site further upwards, as previously 
observed. Both \o an d X±i appear to grow unbounded as u> — > 1 (the small amplitude limit) in 
agreement with the observation that the length of the continuation neighbourhood vanishes in this 
limit. The graph suggests that as to — > (T — > oo), xo vanishes, while x±i tends to a nonzero con- 
stant. Of course, this limit is numerically inaccessible, so these remarks are necessarily speculative. 
One can estimate Xn(0) from the numerical work of section ^ by comparing with & a after the 
first continuation step. The results of this calculation are also presented in figure 4 (broken curve) 
for purposes of comparison with the linear method. 

6 Conclusion 

In this paper we have adapted the methods of MacKay and Aubry to prove existence of breathers 
in the TDSG system, despite the fact that there is no limit in which the system decouples into 
independent oscillators (so the existence theorems so far formulated [jl], || do not directly apply). 
Numerical analysis has revealed that the persistence of these breathers away from the h = 2 limit 
depends critically on their period. One reason for this is that the coupling at which the breather 
enters into resonance with the phonons depends on its frequency. However, for most frequencies the 
continuation actually fails before it reaches the phonon resonance region. The underlying mechanism 
for these early bifurcations is at present unknown. A full bifurcation analysis may prove informative. 

In no case does existence persist for all h £ (0,2], so the breathers found here are unconnected 
with the breathers of the continuum sine-Gordon system (which the TDSG system approaches as 
h — ► 0). Indeed, they look qualitatively different from continuum breathers in that they have 
sombrero rather than hump shaped initial profiles. In this respect they appear to differ also from 
breathers found in previous studies of spatially discrete systems. By means of a simple linear method, 
we have confirmed the uniform trend towards sombrero-shaped breathers in the TDSG system. It 



10 



would be interesting to apply the same linear method to a variety of more conventional oscillator 
networks, since this allows one to examine the direction of continuation of one-site breathers away 
from the weak coupling limit with a minimum of computational effort. 
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Appendix 



The proof of invertibility of DFq rested on the following lemma, whose proof we now give. 

Lemma |3] For all a G (0, 1), all nontrivial solutions of the linear difference equations 

An 2 — An + a x n — 2 Al 
a(2n-l)(n + l) " + n-1 



/N.i .i .i 'Hi ~ till ~r u ,i ii — £ . 1 

(a) A\+A\ = 0, A] l+l + - 7 __ T7 __ T 4 1 l + ;— =0 n > 2, 



8n(n — 1) , An 2 — lOn + 5 .., 
(b) Al + ,Al = 0, Al +1 + a(4n2 l +27l V 2 n + 4^ + 2,-1 ^ = ° " ' 2 



2(4n 2 - 1) 4n 2 -6n + l, . 

c QA\ + llaA\ = 0, A 3 n+1 + —\ — - — '—A 3 n + —= — =0 n > 2 

w 1 2 a(4n 2 + 6n+l) ™ An 2 + 6n + 1 

are exponentially divergent. 



Proof: Each of the difference equations (a), (b), (c) is of a type subject to the following theorem of 
Poincare |l7| : 

Theorem 13 (Poincare) Let f±, f% be functions with limits p\,P2 respectively as n — > oo, such that 
the polynomial P(X) = A 2 + piX + P2 has two distinct real roots X%, A2. Then any solution of the 
linear second order difference equation y n +i + fi(n)y n +f2(ri)y n -i = has the asymptotic behaviour 
linin^oo = A, for either i = 1 or i = 2. 

In each of the three cases of interest here, P is the same polynomial, namely P(A) = A 2 + -A + 1. 
For a 6 (0, 1), this has two distinct negative roots whose product is unity. Let Ai denote the larger 
root, so that |Ai| < 1 while | A2 1 > 1. To prove exponential divergence of A 17 , a — 1,2,3, it is 
convenient to define the quotient sequences 

K - (58) 

which satisfy first order nonlinear recurrence relations. Note that the sequence B" exists if and 
only if it never vanishes. We will assume that these sequences exist, and justify the assumption a 
posteriori. By Poincare's Theorem, for each a the sequence P>° — > |Aj| for either i = 1 or % = 2. To 
prove exponential divergence of A a therefore, it suffices to show that B a cannot converge to a limit 
less than 1. 

Part (a): From the difference equation (a), one finds that B 1 satisfies the first order recurrence 
relation 

1 An 2 - An + a n — 2 1 
Bn = a(2n-l)(n + l) "n + lB^ (59) 
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with initial condition B\ = 1. We claim that B\ > 1 for all n £ N. The proof is by induction. 
Assume Bi, > 1. Then 



Bl = 



An 2 -An 1 



a(2n-l)(n + l) (2n-l)(n + l) fj + 15^ 



An 2 — An + 1 n - 2 , , 

> 7 — - — — — = 1. 60 

(2n-l)(n + l) n + 1 

Since > 1 for all n £ N, the sequence exists (never vanishes) and cannot converge to |Ai| < 1. 
Hence B 1 — ► |A 2 | > 1, and A 1 is exponentially divergent. 

Part (b): From the difference equation (b), one finds that B 2 satisfies the first order recurrence 
relation 

p 2 = 8 n(n-l) An 2 - Wn + 5 1 
" a 4n 2 + 2n - 1 4n 2 + 2n - 1 B 2 n _ x [ ' 

with initial condition B\ — i . Since the solution of ( |6l|) depends on the parameter a we will denote 
it B 2 (a), and further, allow a to take all values in (0,1], denoting B 2 (l) by B. The first step is 
to prove that B exists (is nonvanishing) and does not converge to a limit less than 1. Assume 
~B n -i > n/(n + 1). Then for all n > 2 

— 7i + l 8n(n - 1) An 2 - lOn + 5 n + 1 

" ~ n + 2 > 4n 2 + 2ra - 1 4n 2 + 2n - 1 n 
6n - 10 

(4n 2 + 2n-l)(n + 2)n > (62) 

Explicit evaluation of the first few terms yields that B3 = |j > |. Hence, by induction, for all 
n > 3, B n > (n + l)/(n + 2), so 5 exists and does not converge to a limit less than 1. 

Now consider the sequence B 2 (a) for a < 1. We claim that each term B 2 (a) is a nonincreasing 
function on (0, 1]. To see this, observe that B\{a) = J is manifestly nonincreasing, and assume that 
B 2 _ 1 (a) is nonincreasing for some n > 2. Then B 2 _ 1 (a) > i? 2 _ x (l) = > by the paragraph 

above, so B 2 _ 1 (a) is positive and nonincreasing. From the relation ( |6l|) it immediately follows that 
B 2 (a) is nonincreasing, so the claim is proved. Hence B 2 {a) is strictly positive (and hence exists) for 
all a £ (0, 1). Further limB 2 (a) > limB > 1 > |Ai|, so B 2 — * | A2 1 and A 2 is exponentially divergent. 

Part (c): The argument is almost identical to part (b). The sequence B 3 (a) satisfies 

g3 _ 2 4n 2 - 1 4n 2 - 6n + 1 1 

" _ a An 2 + 6n + 1 ~ 4n 2 + 6n + 1 ^ ^ 

with initial condition Bf(a) = jj-. Once again we allow a to take all values in (0, 1] and denote 
B 3 (l) by £?. We claim that for all n £N, B n > n/(n + 1). Clearly, the claim holds true for n = 1. 
Assume B„_i > (n — l)/n for some n > 2. Then 

— 4n 2 — 1 4n 2 — 6n + 1 n 
^ ^ 2 ' 

™ 4n 2 + 6n + 1 4n 2 + 6n + 1 n - 1 

n 2 

^ B " ~ ^TT > (4n 2 + 6n+l)(n+l)(n-l) > °' (64) 

so that B n > n/(n + 1) and the claim is proved. Hence B is strictly positive (hence exists) and 
cannot converge to a limit less than 1. 

We next claim that for each n £ N, B^(a) is a decreasing function on (0, 1]. Clearly B\(a) = 
is decreasing. Assume that B 3 _ 1 (a) is decreasing. Then Bf l _ 1 (a) > B 3 _ 1 (l) = B n -i > 0, so 
Bf l _ 1 is a positive decreasing function of a. It follows from ([33]) that B^ia) is also decreasing, and 
the claim is proved. Hence B 3 (a) is strictly positive (and hence exists) for all a £ (0, 1). Further 
lim£? 3 (a) > limi? > 1 > |Ai|, so B 3 — > |Aa| and A 3 is exponentially divergent. □ 
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